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Abstract: Mixed action theories with chirally symmetric valence fermions exhibit very de- 
sirable features both at the level of the lattice calculations as well as in the construction and 
implementation of the low energy mixed action effective field theory. In this work we show 
that when such a mixed action effective field theory is projected onto the valence sector, 
both the Lagrangian and the extrapolation formulae become universal in form through next 
to leading order, for all variants of discretization methods used for the sea fermions. Our 
conclusion relies on the chiral nature of the valence quarks. The result implies that for all 
sea quark methods which are in the same universality class as QCD, the numerical values of 
the physical coefficients in the various mixed action chiral Lagrangians will be the same up to 
lattice spacing dependent corrections. This allows us to construct a prescription to determine 
the mixed action extrapolation formulae for a large class of hadronic correlation functions 
computed in partially quenched chiral perturbation theory at the one-loop level. For specific 
examples, we apply this prescription to the nucleon twist— 2 matrix elements and the nucleon- 
nucleon system. In addition, we determine the mixed action extrapolation formula for the 
neutron EDM as this provides a nice example of a ^-dependent observable; these observables 
are exceptions to our prescription. 
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1. Introduction 



There has recently been a rapid growth in the use of mixed action or hybrid lattice QCD ^] 
in the numerical computation of hadronic matrix elements []|, ^ || ^, 0, ^, [R^, |ll], [l^, [b|, 

m, n, m, 0, m. in response, there have been significant complementary developments in 

our theoretical understanding of mixed action (MA) lattice QCD through the use of mixed 

,1,0. Mixed 



action effective field theory (EFT) pL 20, P, pL Eg, M 25, 26, 27, 28 



action calculations allow one to use fermion discretization methods which respect chiral sym- 
metry in the valence sector during the construction of the hadronic source and sink operators. 
This is done in the background of numerically cheaper discretization methods in the sea sec- 
tor (which generally explicitly violate chiral symmetry) during the generation of the gauge 
field configurations which contain the dynamical quark-antiquark polarization loops. The 
main motivation stems from the significant numerical cost [j32[ , |33| of simulating either dy- 
namical Kaplan (domain-wall) fermions 



or dynamical overlap fermions 



35 



mm, _ 

in the chiral regime as compared to Wilson fermions |3^] (including clover [4C] and twisted 
mass El]]) or staggered fermions [f42| , E| . Kaplan and overlap fermions are often referred 
to as Ginsparg- Wilson (GW) fermions as they (approximately) satisfy the GW relation 44] 
and thus retain chiral symmetry on the lattice |45f | (modulo the quark masses). The cost of 
these MA calculations employing GW valence fermions is then only the cost of propagator 
generation in the background of the dynamical gauge configurations (which presumably use 
one of the numerically cheaper varieties of fermions for the sea sector). The most popular MA 
scheme |, § §, g, §, g, 0, [jl], 0, 0, [RJ Q, Q was developed by the LHP Collaboration g, g] 
in which domain-wall valence propagators are generated on the asqtad-improved pub- 
licly available MILC configurations pq| . 

In a recent paper pS[| , we showed that the chiral symmetry of the valence fermions 
suppresses sources of chiral symmetry breaking arising from the sea sector such that for many 
mesonic observables, there are no lattice spacing dependent counterterms through next-to- 
leading order (NLO) in the MA EFT. Furthermore, this valence chiral symmetry is strong 
enough to suppress all explicit lattice spacing dependence of these mesonic observables, with 
the exception of modifications to the correlation functions arising from the unphysical hairpin 
contributions. This exception is due to the lack of unitarity inherent in MA or partially 
quenched (PQ) calculations [49, 5C]. These properties are only transparent when one uses an 
on-shell renormalization scheme, expressing correlation functions in terms of lattice-physical 
parameters measured directly in the calculation, such as the pion mass or decay constant, 
or /„-; in much of our discussion below we will assume that this renormalization scheme has 
been utilized. 

In this paper, we build on these results and extend them such that we can formulate 
a convenient prescription for converting quantities computed in partially quenched theories 
to expressions valid in mixed action theories. Our prescription is valid for a wide range of 
observables in a range of mixed action theories. To be clear, we state our requirements here 
and will discuss them more extensively throughout this work. In terms of the mixed action 
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theory, we require that the hairpin structure in the mixed action theory is the same as in the 
partially quenched theory, and we require that the valence quarks are chiral. In terms of the 
observables, we require that there is no dependence on the CP violating 8 parameter. We 



will discuss these requirements below and later in Sec. 3.1 we will describe how the neutron 
EDM, which of course depends on 9 requires a modification of our prescription. 

We begin by observing that mixed action EFTs describing the light mesons have one 
unphysical operator appearing at leading order (LO) which is universal in form, regardless 
of the discretization used in the valence or sea sector, with only the numerical value of the 
coefficient depending on the actions used; the coefficient is known as Cmix m the literature. 1 
In Section |2| we prove that for vertices with 2N mesons from this operator, for which two 
of the mesons are of a mixed valence-sea type, and the rest are purely valence (or sea), this 
operator functions identically to the LO operator involving the quark mass term. Also in 
Sec. ||, we construct the most general MA Lagrangian projected onto the valence sector of the 
theory. This is particularly relevant because for correlation functions computed with valence 
fermions, we only need the counterterms of the valence sector. 

Our results imply that a large cancellation of potential lattice spacing dependent countert- 
erms occurs in MA theories. This means that MA extrapolation formulae have a continuum 
functional form with only slight modifications, when expressed in lattice-physical parameters. 
Combined with our work in Section |2[ this allows us to construct the general prescription we 
alluded to above. The prescription converts PQ chiral extrapolation formulae through the 
one-loop level into the corresponding MA extrapolation formulae, allowing one to make use 
of the extensive literature on partial quenching. As we will discuss in some detail in Sec. 
our prescription requires three key components to be valid; the valence fermions are (approx- 
imately) chirally symmetric modulo the quark masses, the hairpin structure of the theory is 
the same as in partially quenched chiral perturbation theory (PQxPT) |4^, [5(], [5^] and 
the 8 term is negligible. 

In Section 0, we explicitly determine the MA formulae of several observables which are 



non-trivial examples of our prescription and are of current interest; Sec. 3.1, the neutron 



EDM which provides a nice example and requires a slight modification of our prescription, 



in Sec. 3^ nucleon twist-2 matrix elements and in Sec. [3.3| , nucleon-nucleon scattering. In 
Section || we comment on our results and conclude. In the appendix, we describe why mixed 
action theories involving a twisted mass sea satisfy the requirements of our prescription. 

2. MA Effective Field Theory 



Mixed action EFT is a natural generalization of PQxPT J49L 50, 51, M] reducing to it in 



the continuum limit. 2 Partially quenched xPT is constructed from the underlying theory, 



partially quenched QCD (PQQCD), analogously to how xPT |53|, |5J, [55| is constructed from 



1 This parameter, or equivalently the mixed meson mass renormalization has recently been calculated for 
domain- wall valence and the coarse MILC staggered fermions JlTj] . 

2 For an introduction to MA EFT see Refs. |l9| |f, §§§3], |§. 
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QCD. Partially quenched QCD exhibits an approximate graded chiral symmetry (for light 
quark masses, tuq <C Aqcd)> with N v valence and N v ghost quarks and N s sea quarks, 3 

SU(N V + N S \N V ) L ® SU(N V + N S \N V ) R ® U(l) v , 

which is explicitly broken by the quark mass terms, niQ. It is then assumed, as with QCD, that 
the vacuum of PQQCD spontaneously breaks this symmetry down to the vector subgroup, 
giving rise to the PQ pseudo-Goldstone modes. The PQxPT Lagrangian is then constructed 
with a spurion analysis such that all operators respect the symmetries of PQQCD. These 
symmetries are further broken explicitly by mixed action effects. At finite lattice spacing, a, 
there is no symmetry which mixes the valence and sea quarks, breaking the symmetry to a 
direct product of valence and sea sectors, 

SU(N V + N S \N V ) L ® SU(N V + N S \N V ) R 

SU(N V \N V ) L ® SU(N V \N V ) R SU(N S ) L ® SU(N S ) R . (2.1) 

a^O 

For sufficiently small lattice spacing compared to the non-perturbative scale, these effects are 
perturbative and can be treated in an EFT framework. All operators in the MA Lagrangian 
which do not explicitly depend upon the lattice spacing are given by their PQ equivalents, 
with the same value of the corresponding low energy constants (LECs). There will be new 
operators with explicit lattice spacing dependence, some of which arise from the mixed action 
effects and break the PQ symmetry, and whose exact form depends upon the lattice actions 
used. These new unphysical operators will contribute to correlation functions of observable 
quantities, for which the extrapolation formulae can be determined from the appropriate 
mixed action EFT and then used to remove these unphysical contaminations from MA lattice 
QCD calculations. 

However, mixed action theories that have chirally symmetric valence fermions, such as 
domain-wall fermions in the infinite 5 th dimension limit, or overlap fermions with a per- 
fect overlap operator, give rise to chiral extrapolation formulae for a large class of valence 
observables which are identical in form through next-to-leading order, with any and all dis- 
cretization methods used for the sea fermions. Provided the various sea quark discretization 
methods are in the same universality class as QCD, and that the lattice spacing dependent 
chiral symmetry breaking is perturbative 4 , that is, the lattice spacing a is small compared to 
the scale of chiral symmetry breaking A x in the sense that oA x <C 1, in the same way that 
the quark masses m q give rise to perturbative chiral symmetry breaking since m q <C A x , the 
only difference between these various mixed action theories will be in the numerical value 

3 See Ref. Ji2| for a complete discussion of the PQ symmetry group. 

4 Notice that if the lattice spacing dependent sources of chiral symmetry violation are too large, then chiral 
perturbation theory is simply not the correct effective field theory describing the low energy dynamics of the 
underlying lattice theory. Similarly, chiral perturbation theory would not be relevant in nature if quark masses 
were all large compared to the scale of chiral symmetry violation. 
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of the unphysical counterterms which enter the chiral extrapolation formulae. Furthermore, 
these extrapolation formulae are sufficiently continuum like, due to the good chiral properties 
inherited by the chiral symmetry of the valence fermions, that they can be determined from 
the corresponding formulae in partially quenched xPT. In the rest of Sec. |2|, we present the 
formalism necessary to understand this claim and then provide our mixed action prescription. 



2.1 Matching the 0{a 2 ) operators 

To construct the mixed action effective Lagrangian, one must first construct the continuum 
Symanzik quark level Lagrangian [156], 57 1 which respects all the symmetries of the underly- 



ing lattice action. Then one performs a spurion analysis on this continuum Lagrangian to 
determine the operators in the mixed action EFT |58|]. 5 A specified power counting orders 
the infinite tower of operators entering the low energy Lagrangian. In this work we consider 
the general small parameter to be 

2 

e ~ ~a A x . (2.2) 

X 

The order at which the LO lattice spacing dependent operators appear depends upon the spe- 
cific action used. However, to determine the counterterm structure of the chiral Lagrangian 
relevant to valence quantities through NLO in the mixed action EFT, we only need to un- 
derstand the quark level operators of the valence and mixed sectors of the theory. The sea 
quark operators will be important for determining the additive mass renormalization of the 
sea-sea mesons, which is important for understanding the hairpin interactions, but otherwise 
the sea quark Symanzik operators will only lead to trivial renormalizations of counterterms 
relevant to valence quantities. This will become more clear in our discussion below. We 
also stress that this only holds through the leading loop order in the EFT, after which the 
extrapolation formulae for valence quantities will become dependent upon the details of the 
underlying lattice action in the sea sector. 

The chiral symmetry of the valence sector prohibits operators of dimension five in the 



Symanzik action [19, 2C, 22|. Therefore, we begin with the dimension-6 operators which mix 
the valence and sea fermions, and are therefore necessarily four quark operators. There will 
be 0(a?) operators in the valence sector but they will not break chiral symmetry because 
of the good chiral properties of the valence fermions. These will then be indistinguishable 
from the operators already in the chiral Lagrangian and will amount to renormalizations of 
the physical operators. There are also operators which break Lorentz symmetry, but are 
singlets under the hypercubic group which we do not consider here as they are generally 
higher order than we are working to. In Ref. [21], these Lorentz violating operators have 
been analyzed for the baryons. As we mentioned above there will also be 0(a 2 ) operators in 



We note that it is unnecessary to consider the color structure of these effective continuum quark level 
operators as this does not play a role in the pattern of symmetry breaking used to construct the EFT operators. 



- 5 - 



the sea-sector, 6 but these transform as singlets under chiral rotations of the valence fermions, 
and thus can be absorbed into generic 0(a 2 ) operators needed to renormalize purely valence 
correlation functions. 7 Therefore, to construct the mixed action chiral Lagrangian we only 
need to consider the mixed valence-sea operators at this order. Much of this discussion can 
be found in Refs. 



19 



21, 22, 24 



The valence-sea mixing Lagrangian at dimension six for chirally symmetric valence fermions 
and any type of sea fermion is given by four-quark operators which are products of valence 
and sea quark-bilinears that independently respect chiral symmetry in the valence and sea 



sectors respectively |T9], 2C, 21, 22, 24 1. However, because these operators explicitly break 
the partially quenched symmetry relating the sea and valence fermions, despite the fact that 
they are constructed from chirally symmetric quark-bilinears, they give rise to additive mass 
corrections for mixed hadrons composed of both valence and sea fermions. This, for exam- 
ple, is how a mixed pion of domain-wall valence and sea quarks, but with different values of 
the Wilson- "r" parameter, or different 5 th dimensional extent, are subject to additive lattice 
spacing dependent mass corrections. In the non-mixed action limit, these operators are no 
longer allowed by the symmetries of the lattice action and must vanish. At dimension-6, there 
are only two allowed mixed action operators which are both 4-quark operators, 



4L = « 2 ^Mix [Ql^Vv Q) (Q-fVs Q) + a 2 C^ lx (Q W V V Q) {Qf-ysPs Q) 

= « 2 C M f x (Ql l^v Ql) {Ql l»Ps Ql) + a 2 C M f x {Ql IpPv Ql) (Qr l^s Qr) 



R] 



(2.3) 



where Vy and Vs are valence and sea projection operators respectively. The coefficients 

Cyk due to parity conservation. The introduction of Vh,a\ and 



riLL r<RR ot ,j r<LR 

°Mix — U Mix ancL °Mix 

Vyr S \ is convenient for spurion analysis after which one can set "Pyrgs 
Under chiral transformations, Ql 
under these transformations if 



LQl and — ► RQr. Equation (E 



V(S) 

7V(S)- 
will be invariant 



V R 

' V(S) 



V V(S) 



RV v{s) r] - 



(2.4) 



In MA^PT, the hadronic fields transform as in PQxPT under the chiral transformations, 



which are ||, |g |l|, |2| 



S(x) -» LS(x)i? f , i{x) ^ L£{x)U ] {x) = U{x)£{x)R ] 
Bijk (-l) 7?i ^' +7 ' m ) + ^ + ^^' 7fc+7 '")c/(x)/C7(x) i -C/(x) fe n ^ mri: 



(2.5) 



6 For unimproved Wilson sea fermions, there will be 0(a) operators. This leads to an extra complication 
with the hairpin interactions, but otherwise does not modify our prescription. We will address this in more 



detail in Sec. 2.2 



7 Of course there are exceptions to this rule, but these occur only in special cases which are related to the 
class of quantities which do not follow our discussion and prescription. We will postpone the discussion of 



these special cases to Sec. 3.1, where we will discuss them in the context of the neutron EDM. 
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where £ contains the meson fields, £ = ^ is a spin-1/2 baryon field and C/(x) is a 
complicated transformation which depends upon the mesons, and thus on spacetime. The 
chiral transformations of the spin-3/2 fields, T are identical to those of B and the heavy 
meson transformations can be found in Ref. [BH]. The grading factors, rji keep track of the 
(anti) commuting nature of the different quark fields, where the sea and valence quarks are 
anti-commuting and the ghost quarks are commuting. For our discussion we will be mostly 
interested in products of fields with purely valence fermions as these will be the relevant 
degrees of freedom to construct the counterterms for valence quantities at NLO. With this 
restriction, it is often more convenient to think about how the valence fields transform. The 
valence-nucleon field Ny, for example, can be related to the valence projected B field and 
transforms more simply under the reduced chiral transformations 

B ijk = —7= (eijN k + e ik Nj) for i,j, k G valence (2.6) 
V6 

N v -» U(x)N v . (2.7) 



The projection operators transform as in Eq.( |2.4D and thus we also have 
which allows us to define the projectors 



with 

^(SJ-^fl)^. ( 2 - 10 ) 

which are even (+) and odd (-) under parity. Now we can construct the set of 0(a 2 ) operators 
in the MA chiral Lagrangian relevant for valence quantities, which must have the same chiral 
symmetry properties as Eq. ( |2.3| ), 

£Mix = a 2 (v Mix + v Mix + v Mi^ (2 n) 

The mixed action meson potential was first determined in Ref. ||20| ] . We will deduce the 
potential here on the basis of symmetry considerations in order to emphasize its universality. 
The symmetry breaking of Eq. (2.1) enlarges the number of operators which we can form 



from the meson field X, and which are invariant under the symmetry. These operators occur 
at order C(a 2 ), which we take to be of order m q in our power counting. Consequently, the 
lowest order operators which are invariant under the reduced symmetry are 

0i = a 2 str{V s ZVs^), 2 = a 2 st?(V v Ws^), (2.12) 
3 = a 2 strCP 5 £?V£t), q 4 = a 2 str{V v ZVv^). (2.13) 
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However, using the identities Vs + Vy = 1 and X^E = 1, it is easy to see that these opera- 
tors are equal up to constant numeric factors, so that there is only one non-trivial operator 
invariant under the reduced symmetry group at this order. Thus, the mixed action meson 
potential is conventionally given by 

Vf ix = C Mix str (T 3 £T 3 £t) (2.14) 

with 

T 3 = V s - Vy ■ (2.15) 

Of course, the quark level Lagrangian given in Eq. fl2.3| ) is universal by the same reasoning. 
Note, this also applies to unimproved Wilson sea fermions as well. 

There are additional 0(a 2 ) mixed contributions involving the sea quarks. These operators 
are singlets under chiral rotations of the valence quarks so they can be absorbed into the PQ 
coefficients and, in general, will give additive a 2 contributions to the PQ coefficients. For 
example, 

V^ ix = C' N {NyV+Ny) str[P+] 

— ► C N N V N V , (2.16) 

which we see is an additive correction the valence nucleon mass. 8 Therefore, we can just 
focus on the 0(a 2 ) contributions of the valence fields from now on as the contributions from 
the mixed terms are indistinguishable under valence-chiral rotations. If we are interested in 
matrix elements with no pions in the external states, the nucleon counterterms behave as if 
C = C 1 an d therefore 

^W^Vw. V ns)^^ ( 2 - 17 ) 

immediately eliminating all operators with insertions of 'Py(s) f rom our consideration. Sim- 
ilarly, operators of the form NyVyVg Ny also vanish by projection. The above discussion 
generalizes to any single nucleon 0(a 2 ) counterterm of the form 

NyOVvN v sti[Vs] , (2.18) 

and therefore we can parameterize the single nucleon counterterms relevant at NLO with the 
one operator given in Eq. ( |2.16 ). It is important to note that this term, str [Ps], respects 



chiral symmetry. 

The mixed action two-nucleon (two-baryon, two-heavy meson) Lagrangian can be con- 
structed in a similar fashion. The mixed two-nucleon potential contains only two operators 
relevant at NLO, 



Vnn = C a NN (NyNy) + C b NN (N v S»N v y 

t S0) (NTpf^Ny)\NTpf^Ny 



Di Sl) (N^Ny) 1 (N^Ny) , (2.19) 



3 The operator (NvPg Nv) strp 3 ^] does not contribute because str["P^] = 0. 
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where P^ are the S- wave projection operators for channel-j in the two nucleon system [63]. 
Through the order we are working, there are no counterterms for the other angular momentum 
projections because only the S'-wave two- nucleon wavefunctions are non- vanishing at the 
origin, and so counterterms for the higher partial waves must also contain derivatives pushing 
the mixed action counterterms beyond the order we are working. 

Prom the above construction, one observes that the leading 0(a 2 ) effects in V^f lx and 

c 

NN 



V^fe appear to additively renormalize the sea quark mass dependent terms 



N v N v str[m Q ] , (2.20) 
D'f" ( N$pv>N v V (n^P^Nv) str[S m Q ] , (2.21) 



and 

VO') ( atToU) at.,V ( ktT.pU) 



in the PQxPT Lagrangians [61, |62[ and [64] respectively. However, it is important to stress 
that these lattice spacing dependent operators do not break chiral symmetry, and will differ 
from the sea quark mass dependent operators at higher orders. These terms are allowed 
because the nucleon mass and the NN interactions are not protected by chiral symmetry. 
As discussed in Ref. [29], the 0(a 2 ) correction in V^ lx only gives additive renormalization 



to the valence-sea meson masses while the subleading 0(a 2 m q ) corrections only give additive 
renormalization to the LO meson decay constant / and the chiral condensate, Bq. Thus 
one can obtain the mixed action EFT results at 0(a 2 ) from the PQxPT results almost 
effortlessly. 9 We re-emphasize that this entire discussion and the prescription we present 



in Sec. 2.5 is dependent upon these very benign mixed action lattice artifacts. This simple 
behavior does not hold beyond the leading loop order except in special cases. With these 
caveats in mind, we now provide the mixed action Lagrangian relevant for determining the 
chiral extrapolation formulae for all mixed action theories with chirally symmetric valence 
fermions. 

2.2 The Mixed Action Lagrangian 

For no reason other than author bias, we present the Lagrangian in Minkowski space-time, 
despite the lattice theories being constructed in Euclidean space-time. The LO mixed meson 
Lagrangian is 10 

C T A) = Y str ( 5 M S ^ st ) + 4r Stl " ( mQSt + Sm k) + a2CMixStr ( T 3 ST 3S f ) + a 2 V sea . 

(2.22) 

We have assumed that the sea-quark action is either 0{a) improved, or has scaling violations 
(lattice spacing artifacts) beginning at 0(a 2 ) or 0(a s a 2 ). This is not essential for our discus- 
sion, but as we will briefly discuss below, it simplifies the structure of the hairpin propagators 

9 This entire discussion can also be extended to heavy mesons as well as baryons with a heavy quark. These 
fields also transform with U (x) under chiral transformations and so the construction of the Lagrangian is very 
similar to that of the nucleon presented here. The relevant partially quenched Lagrangians can be found in 
Ref. [^o| for the heavy mesons and in Refs. J65| , prj] , |67| for the baryons with heavy quarks. 
10 See Refs. |2q, M for our conventions labeling the meson fields and quark masses. 
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such that to the order we are working, they are identical to those of PQxPT |49| , 5C, 51, 52]. 
We will nevertheless adopt this assumption below. 

The LO Lagrangian, Eq. ( |2.22 ), gives rise to the masses of the various pseudo-Goldstone 
mesons, with the LO masses for a meson composed of valence quarks, v, sea quarks s or an 
admixture given by 



m 



viv 2 

™2 



B (m Vl + m 



"V%) 1 



m vs = B (m v + m s ) + a A M i x 
m 2 = B (m Sl + m S2 ) + a 2 A sec 



(2.23) 



with 



a 2 A 



Mix 



. 16Cmi 
~1 2 ~ 



(2.24) 



and a 2 A sea determined from a 2 V S ea- The LO Lagrangian also leads to the well known double- 



poles or hairpin propagators amongst the flavor diagonal mesons [49, 50, 51, 52]. The mo- 
mentum space propagator between two flavor diagonal mesons of quark type a and b is given 
by 



ab 



ih z k + ie) 



P 



mi 



+ ie 



N 



f (p 2 - m 2 a + ie)(p 2 - m 2 b + ie) Uk'LiiP 



where 



+1 for a = valence or sea quarks 
— 1 for a = ghost quarks . 



ml, + ie) 



(2.25) 



(2.26) 



In Eq. ( 2.25 ), k runs over the flavor neutral states {4>jj, • • • , 4>rr) and k' runs over the mass 
eigenstates of the sea sector, including the additive lattice spacing mass corrections, Eq. ( p. 23 ). 
To help quantify the unitarity violating corrections arising from these double pole propagators, 
we have introduced partial quenching parameters [29], which are differences between the pole 
masses of the sea-sea and valence- valence mesons, 



A 2 



™ 2 jj 



m: 
m 



2 

uu 

2 

ss 



■■ 2B (mj - m u ) + a 2 A sea + . . . 
2B (m r - m s ) + a 2 A sea + . .. 



(2.27) 



where the dots denote higher order corrections to the meson masses. Using these PQ pa- 
rameters, one can rewrite the hairpin propagators in a particularly simple form allowing for 
a transparent identification of the unphysical unitarity violating contributions to correlation 



functions arising from the hairpin interactions [25, 29]. 

For discretization errors in the sea sector which begin at O(a), such as Wilson fermions 
there are two modifications we need to make to A 2 U and A 2 S . First, the lattice spacing cor- 
rections to the sea-sea meson mass begin at O(a). Second, there are additional hairpin 



interactions whose coefficients depend upon the lattice spacing [23]. This is not problematic 
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to our prescription because these extra hairpins can be treated as an additional additive 0(a 2 ) 
correction to our partial quenching parameters [E3], 



^ 2 ju — ► ™% ~ m lu + a2 lss N s 



A 2 = A 2 



A 2 = A 2 

i "rs "rs 

a=0 



= 2B (mj - m u ) + aW + a j ss N s 
A 2 S — ► 2B (m r - m s ) + aW + a 2 7ss iV s . (2.28) 

If we work consistently to 0(a 2 ), we must include these extra hairpin effects, even though 
they are formally subleading to the 0(a) term in the partial quenching parameters. However, 
determining 7 SS is difficult because it is an additive mass correction to the rf mass. This shift 
in partial quenching parameters will also invariably cause a shift in the numerical value of un- 
physical lattice spacing dependent counterterms but this is accommodated without changing 
the structure of the extrapolation formulae. We will generally assume that the sea fermion 
scaling violations begin at 0(a 2 ) or higher as this is most relevant for lattice calculations of 
the present and future. We will also commonly use the partial quenching parameters in the 
continuum limit to denote strictly differences in the sea and valence quark masses, 

(2.29) 

o=U 

Mixed Action Single Baryon Lagrangian 

The mixed action Lagrangian for the single nucleon (baryon) and interactions with the pions 
(mesons) is given by 

C% A) = i(Bv DB) + 2a (PQ) (BBM+) + 2^ Q) (BM+B) + 2a^ Q) (BB) str(M+) 
- (T» [iv ■ D — A] T M ) + 2 7 Sf Q) (f " M + T^)- 2a™ (T%) str (M+) 
+ 2a {PQ) (BS^BAfj) + 2(5 {PQ) (BS ■ AB) + 2H (f u S^A^ %) 

+ Jlc (f " ApB + BAp T") + a 2 C^ {BB) - a 2 C T a (f^) . (2.30) 

Notice that there are only two new operators in the mixed action Lagrangian as compared 
to the partially quenched Lagrangian for either SU(6\3) |H| or SU(A\2) |6^]. If we were not 
projecting onto the valence sector of the theory or we were interested in working to higher 
orders, we would need additional lattice spacing dependent operators. For example, there 
are i3- field operators similar to those with the mass spurion field M + with chiral symmetry 
breaking a 2 spurions instead. However, for extrapolation formulae of valence quantities, these 
operators all collapse into the form given in Eq. ( fOol) at the order we are working. 

The flavor structure and contractions of these fields, defined with the braces, ( ), can 



be found in Ref. |61[] for the three flavor EFT and in Ref. [S2| for the two flavor EFT. One 
often encounters formulae expressed with the more familiar %PT couplings instead of the PQ 
couplings, for example the SU{A\2) couplings can be related to the SU{2) couplings, 

9 A = \*t PQ) ~ \^ PQ \ 91 = \*™ + \^ PQ) 

U = ffAA, C = 9 AN (2.31) 
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The mass spurion field is given by 

M+ 

and the axial-meson field is 



2 (£rn Q £ + im ] Q i 



(2.32) 
(2.33) 



Mixed Action Two Nucleon Lagrangian 

Following the normalization and conventions of Ref. 
grangian also involves simple modifications to the continuum two-nucleon Lagrangian 68 



], the mixed action two-nucleon La- 



70, 71, 72, Um and is given by 



t t p(j) ■ 



TT jdU) 



a 



TpU) *r V ( ATTpti) 



D 2A (n^P^Nv^ (N$P t if) 2Bm Q N v 



T T T>(j) 



N^P^Ny) iV^ u; (Vr^v 



,2 n 0') 



--Mix 



(2.34) 

where there is an implicit sum over j = j 1 ^, 3 S'i}. There are additionally operators with the 
j = 3 Di projectors as well as the coupled 3 S±- 3 Di system necessary for understanding the 
deuteron, but as we discussed in Sec. [O] the lattice spacing counterterms for these channels 
are suppressed beyond the order we are considering. 

In all of these MA Lagrangians, we see there are only very benign lattice spacing de- 
pendent operators. This is crucial to the universal nature of the mixed action extrapolation 
formulae and a key component that allows us to construct our universal prescription. Even 
though the detailed form of the Lagrangian at NLO depends upon the type of sea fermions 
employed, when we project onto the valence sector of the theory, all that is necessary to 
discuss the renormalization of valence quantities, we see this remarkable simplification of 
the relevant Lagrangian, as has been discussed in this section. It is also clear that the sim- 
ple form this Lagrangian takes depends crucially upon the chiral symmetry of the valence 
fermions, without which there would be several additional operators as is seen with the Wil- 
|H, |7|, f|| E01 H twisted mass ^ 0, [[|, [7| and staggered [|(], H, H, H chiral 



son 



Lagrangians. We postpone a discussion of the class of observables we know to not follow our 
prescription until Sec. |3.1| , where we use the neutron EDM as an example to highlight both 
the reasons this discussion fails to accurately describe this quantity as well as how to modify 



the prescription we provide in Sec. to account for this class of observables which is more 
sensitive to the sea-sector. 

2.3 Mixed Action Meson Operator as a Mixed Meson Mass 

To complete our discussion on MA Lagrangians at the one loop level, we prove that the mixed 
action meson operator functions exactly as the LO quark mass operator for a process with 
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two valence-sea mesons, and 2N — 2 valence- valence mesons. Let us consider the mixed lattice 
potential of Eq. (gjp . Making use of Eq. flDjD and the equality EE* = 1, one can show 



that up to a constant this is equal to 

a 2 C Mix str (T 3 £T 3 £t) = 4a 2 C Mix str (v s ZV s ^) 



- 4 ° CMix ^( v /J ^ n!(2iV-n)! " ( 2 " 35 ) 

AT=0 v J 7 n=0 v ' 

The only terms in the sum which contribute to vertices with two mixed valence-sea mesons 
come from either n = or n = 2N, for which one can show 

_ 2N /<n„A2N 



C Mix str (T 3 £T 3 £t) _> 8 a 2 C Mix ]T (j ) str (^)f) 

= 4 a 2 C Mix str [p s (s + £t) ] ; ( 2 . 3 6) 

and therefore the ratio of this mixed meson operator restricted to vertices involving two 
valence-sea mesons and 2N — 2 valence- valence 11 mesons to the valence-sea quark mass con- 
tribution of the same vertices is 

a 2 C M i x str [T 3 £T 3 £t] 16a 2 C M ix 2a (n 

~ -p = a ^Mix • {^-Jt) 



(P/4) S tx[B m Q (E + Et)} 



(2N-2)<j, vv 



This is evident in previous determinations of mixed action extrapolation formulae [25, 26, 27 



3Cj. Thus, the effects of this operator at one-loop act simply as a shift in the mixed 



valence-sea meson masses in all vertices and propagators. 
2.4 Mixed action EFT at one loop 

The symmetry structure of the underlying mixed lattice action determines the operators in 
the mixed action chiral Lagrangian. The symmetries enjoyed by the valence fermions are 
different from those enjoyed by the sea fermions. In the class of mixed action theories we 
consider here, the valence fermions only break chiral symmetry through the explicit quark 
mass term. Therefore, in the meson Lagrangian at NLO, the purely valence spurions are 
identical to the spurions in continuum, unquenched chiral perturbation theory, and so the 
valence-valence sector of the NLO mixed action chiral Lagrangian is the Gasser-Leutwyler 
Lagrangian. This is not the case for baryons which have LO lattice spacing operators as we 
have seen in Eq. flT|g) . The sea sector is different. At finite lattice spacing, the sea sector has 



enhanced sources of chiral symmetry violation — for example, there are additional spurions 
associated with taste violation if the sea quarks are staggered, or in the case of a Wilson sea, 
the Wilson term violates chiral symmetry. Consequently, there are additional spurions in the 
sea sector. Of course, these spurions must involve the sea quarks and must vanish when the 
sea quark fields vanish. 



This also holds for vertices with 2N — 2 sea-sea mesons and two mixed valence-sea mesons. 
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In this paper, we work consistently to NLO in the mixed action xPT power counting which 
we have defined in Eq. (|2.2| ). For meson observables, the NLO operators in the Lagrangian 
are only used as counterterms; that is, at NLO one only computes at tree level with the NLO 
operators. Since the in/out states used in lattice simulations involve purely valence quarks, we 
can project the NLO operators onto the valence quark sector of the theory. Consequently, all 
of the spurions which involve the sea quark fields vanish. Since the remaining spurions involve 
the valence quarks alone, we only encounter the symmetry structure of the valence quarks 
as far as the NLO operators are concerned. These spurions only depend on quark masses 
and the quark condensate itself, and so there can be no dependence on lattice discretization 
effects arising in this way. The exception to this argument arises in the case of double 
trace operators in the NLO chiral Lagrangian; in these cases the valence and sea sectors 
interact in a flavor-disconnected manner, unlike the mixed operator in Eq. ( [2.14 ). If one 
trace involves a valence-valence spurion while the other involves a sea-sea spurion, then the 
trace over the sea may still contribute to a physical quantity, for example the meson masses 
and decay constants. Note that the valence-valence operators which occur in these double 
trace operators must be proportional to one of the two operators present in the LO chiral 
Lagrangian, Eq. fl2j2| ). Thus, as it was argued and demonstrated in Ref. |p9| , for meson 
scattering processes, the dependence upon the sea quarks from these double trace operators 
can only involve a renormalization of the leading order quantities / and Bq. Both the explicit 
sea quark mass dependence and the explicit lattice spacing dependence are removed from 
the scattering processes expressed in terms of the lattice-physical parameters since they are 
eliminated in favor of the decay constants and meson masses which can simply be measured on 
the lattice. When expressed in lattice-physical parameters, there can be no dependence upon 
the sea quark masses leading to unphysical PQ counterterms and similarly there can be no 
dependence upon an unphysical lattice-spacing counterterm. This argument generalizes to all 
meson quantities which are protected by chiral symmetry. For completeness, we summarize 
the discussion in Ref. |29| here. 

The NLO Lagrangian for mesons describing the valence and sea quark mass dependence 
is the Gasser-Leutwyler Lagrangian with traces replaced by supertraces: 



£-GL = Li 



str ( dulld^)} + L 2 str ( d u Zd^) str ( d^d"^ 



+ L 3 str ( djld^d^d^ ) + 2B L 4 str ( dj^d^ ) str 



+ 2B Q L 5 str 



+ Em! 



+ ±B> L 6 



str ( m q T^ + 'Erriq 



+ 4B 2 L 7 



str ( mjs 



m. 



+ 4i?n L$ str ( m q T^m a T) + £m£ SmJ 



(2.38) 



Having a concrete expression for the Lagrangian, we can show explicitly how the sea quark 
mass dependence disappears. The key is that when constructing NLO correlation functions 
of purely valence quarks, we can replace the mesonic matrix $ in the NLO Lagrangian by a 
projected matrix 

$ ^ Vy^Vy (2.39) 
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where Vy is the projector onto the valence subspace. The only sea quark mass dependence 
comes from two operators 

5C GL = 4P L 4 str (d^P v d^P v ^j stv(m q ) 

+ 16B 2 L 6 str (m q Y)Py + P v Y>m) q ^j str(m ? ) , (2.40) 
which leads to a renormalization of the LECs / and Bq 

f^f + 32L 4 B str(m,), f 2 B / 2 P + 64L 6 P 2 str(ro,). (2.41) 

Since the parameters / and Bq are eliminated in lattice-physical parameters in favor of the 
measured decay constants and meson masses, we can remove the dependence on the sea quark 
masses by working in lattice-physical parameters. Analogously, we can remove all the explicit 
lattice spacing dependence. The general MA Lagrangian involving valence-valence external 
states at 0{e 2 n e 2 a ) can be reduced to the following form 

5C M A = a 2 L™ 2 q str (m q P v ^P v + P v ^P V m q ) str (g(P s XP s ) J(P s tfP s j) 

+ a 2 L 9 a2 str (d^P v d^P v ) str (/(P S £P S ) f'(Ps^P s )) + h.c, (2.42) 

where the /'s and g's are functions dependent upon the sea-quark lattice action. These then 
lead to renormalizations of the LO constants, 

f ^ f + 8a 2 L d 2 str ( /( p 5 ip 5 ) f'(P S lP S )) , 

/ 2 P ^ / 2 Po + 4a 2 L™ 9 str ( 5 (P 5 1P 5 ) 5 / (P S 1P 5 )) , (2.43) 

and just as with the sea quark mass dependence, expressing physical quantities in terms of 
the lattice-physical parameters removes any explicit dependence upon the lattice spacing. 

Together, these results show that at NLO, the only counterterms entering into the extrap- 
olation formulae for mesonic observables protected by chiral symmetry are the same as the 
continuum Gasser-Leutwyler counterterms entering at NLO. This lack of unphysical coun- 
terterms is desirable from the point of view of chiral extrapolations, but it also has another 
consequence. Loop graphs in quantum field theories are frequently divergent; there must be 
a counterterm to absorb these divergences in a consistent field theory. Since there is no coun- 
terterm proportional to a 2 or the sea quark masses, loop graphs involving these quantities are 
constrained so that they have no divergence proportional to a 2 or the sea quark masses. This 
further reduces the possible sources of sea quark or lattice spacing dependence. For example, 
mixed valence-sea meson masses have lattice spacing shifts, so there can be no divergence 
involving the valence-sea meson masses. 

2.4.1 Dependence upon sea quarks 

At NLO in the effective field theory expansion, mesons composed of one or two sea quarks 
only arise in loop graphs. In particular, the valence-sea mesons can propagate between 
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vertices where they interact with valence-valence mesons; these interactions involve the LO 
chiral Lagrangian ( |2.22| ). Because the mixing term, Eq. (2.14) is universal, these interaction 



vertices are the same for all discretization schemes provided LO chiral perturbation theory 
is applicable. The sea-sea mesons only arise at NLO in hairpins. Therefore, we see that our 
NLO extrapolation formulae only depend on the LO chiral Lagrangian to quadratic order in 
the sea-sea sector and the LO chiral Lagrangian (with the mixing term) in the valence-sea 
sector. The mixed meson splitting has recently been computed for domain-wall fermions on 



the coarse MILC lattices [17] and was found to be 

a 2 A Mix ^ (314 ± 4 MeV) 2 (2.44) 



for a = 0.125 fm |17j. It has also been determined on the fine MILC lattices as well [p4| . 

Note that the impact of using different sea quark discretizations in our work comes only 
from the value of the LECs of the unphysical operators in the Lagrangian. Therefore, the 
same NLO extrapolation formulae can be used to describe simulations with different sea quark 
discretizations, provided that the appropriate mass shifts are taken into account. In the case 
of staggered sea quarks, the sea-sea mass splitting which occurs in the MA formulae is that of 
the taste-identity, which has been computed [J85|| , and for the coarse MILC lattices, is given 
by 

a 2 A sm = a 2 A 7 ~ (450 MeV) 2 , (2.45) 

for a ~ 0.125 fm. These mass shifts can only appear through the hairpin interactions at 
this order. These terms will generally be associated with unphysical MA/PQ effects which 
give rise to the enhanced chiral logarithms as well as additional finite analytic dependence 
upon the sea-sea as well as valence- valence meson masses (and their associated lattice spacing 
dependent mass corrections). 

For heavy baryon and heavy meson observables, the leading loop corrections are typically 
non-analytic in the quark mass, for example the mass corrections are C(e 3 ). Therefore, the 
only counterterms needed to renormalize the NLO corrections to heavy baryon and heavy 
meson quantities will be those which appear in the LO Lagrangian, Eq. fl2.30| ). The NLO 
loop contributions for these observables will typically involve all types of mesons, valence- 
valence, valence-sea and sea-sea. Therefore, working to NLO, we only need to know the mass 
corrections to the valence-sea and sea-sea mesons. 



2.5 Prescription for Mixed Action Extrapolation Formulae 

We now have all the ingredients to construct a prescription to determine all MA extrapolation 
formulae at the one loop level from the corresponding PQ^PT expressions, for mixed actions 
with chirally symmetric valence fermions and any type of sea fermion. This prescription is 
more useful if one expresses the extrapolation formulae with on-shell renormalization. We 
also restate that this prescription is relevant for theories with chirally symmetric valence 
fermions and the hairpin structure of PQxPT. Given a PQ extrapolation formula make the 
replacements: 
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1. meson and quark masses: exchange the one- loop valence- valence meson masses with the 
lattice-physical meson masses, m uu — ► m^, where m n is the lattice-physical pion mass (or 
appropriate meson mass for other valence quark flavors). Replace tree- level meson masses 
(equivalently quark masses) with the lattice-physical pion mass at the appropriate value of 
the quark mass, with NLO adjustments as needed for consistency in the chiral expansion; for 
example 2B m u -> m 2 - 5m 2 (NLO); 2B m s -> 2m\ - m 2 - 25m 2 K (NLO) + 6ml(NLO). 

2. decay constants: for the LO decay constant / — * /„- — 5f n (NLO) where /„■ is the lattice- 
physical pion decay constant measured on the lattice and 5f 7T (NLO) is the one loop correction 
to this LO value which is entirely determined in terms of the lattice physical parameters. Obvi- 
ously for expressions which are already expressed in the on-shell renormalization, /nlo fn- 
Equivalently, use fx or some linear combination of f n and fx with appropriate NLO adjust- 
ments. 12 

3. mixed mesons: m 2 u — > m 2 u = \m 2 uu + |m|,- + a 2 AMi X5 for a meson composed of a valence 
and sea quarks u and j. 13 

4- sea-sea mesons: m 2 r — > rh 2 r = m 2 r + a 2 A sea for a sea-sea meson composed of sea quark 
flavors j and r with the appropriate additive mass renormalization for a given sea quark 
discretization method. 1 

5. lattice spacing dependent counterterms/ higher dimensional operators: Add lattice spacing 
dependent counterterms (higher dimensional operators) when necessary. Often, this can be 
determined by enforcing the renormalization-scale independence of a given observable. 

3. MA Extrapolation Formulae 

To demonstrate the ease with which our prescription can be applied to the existing partially 
quenched literature, we determine the mixed action extrapolation formulae for several physical 
quantities which are currently of significant interest to the physics community. We note that 
there are several mixed action EFT papers already in existence which provide further non- 



trivial examples of this prescription |2J, |2^, |27|, |28|, |29|, |31|. 

12 This replacement also holds for all couplings which appear in a given formula. For example, LHPC has 
determined the coupling <?a which appears in the nucleon-pion Lagrangian. For an extrapolation of the nucleon 
mass with the same lattice action, one should use this value of qa in the extrapolation formula. 

13 There is no unique way to define the mixed meson mass renormalization; however, the different methods 
only differ at NLO and higher in the mixed meson mass, and therefore this difference will be NNLO or higher 
for all other quantities. As an alternative, for the mixed "kaon" mass, one could make the replacement, 
rn-ru = | m ls" + | m jr + a 2 AMix at the degenerate sea-valence quark mass point. 

14 For a Wilson sea, the mass correction will be linear in the lattice spacing, 8m 2 — aW. For a clover- 
improved Wilson sea ^o), the mass correction is quadratic in the lattice spacing and is given in Ref. po| . 
For a twisted mass sea at maximal twist j4l], it is the tt^ meson mass which enters this expression and the 
mass correction can be found in Ref. jftl], see the Appendix for details. For a staggered sea it is the taste- 
identity meson mass which enters at this expression, which has been measured on the coarse MILC lattices Js5| , 
a 2 A! ~ (446 MeV) 2 . 
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3.1 Neutron Electric Dipole Moment 

The neutron electric dipole moment (EDM) is of great interest both theoretically as well as 
experimentally. A non-vanishing neutron EDM would be direct evidence of CP violations 
which could stem from the #-term in the QCD Lagrangian. There have been several lattice 



calculations of the neutron EDM over the years [86, 87, 88, 89, |90| with continually improving 
techniques and precision. All the calculations to date and the foreseeable future require 
extrapolations to the physical point and given the strong possibility of a mixed action lattice 
calculation of this quantity, it is very relevant to determine the mixed action extrapolation 
formula for the neutron EDM. 

Furthermore, the neutron EDM is interesting because it is an example of a quantity which 
does not follow our prescription. The reason for this is straightforward; the neutron EDM 
is directly proportional to the QCD #-term and therefore is a quantity which is sensitive to 
the axial U(l) chiral anomaly which so to speak "lives" in the sea-sector. This is simple to 
understand in PQ^PT. Upon performing a chiral U(1)a rotation on the valence fermions, 
one must perform an equal rotation upon the ghost "fermions" . In this manner, the change 
in the measure of the anti-commuting valence fields is exactly cancelled by a change in the 
measure of the commuting ghost fields, leaving the theory invariant. 15 However, a chiral 
U(1)a rotation of the sea-fermions is connected to the desired #-term of QCD, hence the 
abuse of language, "the chiral-anomaly lives in the sea-sector." 

This has non-trivial consequences upon the structure of the extrapolation formula for 
the neutron EDM. We can conclude that in the continuum limit, the neutron EDM must be 
proportional to the sea-quark masses (strictly speaking a product of the sea-quark masses), 
because if one of the quark masses were zero, the 0-term is non-physical and thus the EDM 
must vanish in this limit. This rules out counterterms to the neutron EDM which are pro- 
portional to only the valence quark masses. Away from the continuum limit, things are more 
involved and in fact we will need additional operators which we did not include in Eq. ( p.3Cj ). 
We begin with the QCD Lagrangian including the 0-term, 

C = q [zip- m q ] q - \f, u F^ + • (3-1) 

In the continuum limit, the theta-term can be rotated into the quark mass matrix with an 
axial U(l) transformation and then mapped into the chiral Lagrangian. However, at finite 
lattice spacing, this U(1)a transformation will also modify the irrelevant operators in the 
Symanzik action which break chiral symmetry, for example the chromo-magnetic term in the 
Symanzik Wilson Lagrangian will pick up a phase, 

qa^F^q — ► q L e^a^F^q R + q R e* >o ^ 'q L , (3.2) 

with the flavor matrix (j) = diag(^ u , (f)^, . . . ), similar to the quark mass matrix, and 9 = — 
^2j4>j- This will then give rise to lattice spacing dependent operators in the nucleon La- 
grangian which contain this complex phase. For example, the heavy baryon Lagrangian for 



15 This is equivalent to the discussion in Ref. 
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the nucleoli fields will have an operator [|75| , 

C e w D 2a a NaW e + N , (3.3) 

with 

flW j = £^f^ + ^t) . (3.4) 



2 

The aW+ field is parity even and therefore in the absence of the #-term, only contains even 
numbers of pions. However, with the complex phase present, this spurion field also contains 
pions of an odd number, and in particular can contribute to the neutron EDM in the one-loop 



graphs displayed for example in Fig. 1 of Ref. [ 92 ] , in place of the quark mass spurion of the 



nucleon Lagrangian. With mixed action theories, the four quark operators of Eq. ( [2.3D do not 
break chiral symmetry and are therefore invariant under the U{1)a transformation. Therefore 
these mixed operators do not contribute to the neutron EDM. 

In our construction of the mixed action chiral Lagrangian, we have not included certain 
operators which are important in the study of the neutron EDM or any other quantity sensitive 
to the chiral anomaly, which do not generally contribute to observables at this order. They 
stem from four-quark operators constructed from sea-quarks and in this case which also break 
chiral symmetry. To understand these operators, it is convenient to construct chiral symmetry 
breaking spurions of definite parity 

V x,'i = \ (£ f ?W¥ ± iVse-^i) . (3.5) 

In particular, for an 0(a) improved sea fermion action, there are two additional operators 
we should add to the Lagrangian which are important at this order in the presence of the 
0-term, 16 

4T' - 47' 51 + {*r$B) + ^ (bbv£!) Mri!) ■ (3.6) 

We can immediately understand why these operators are not important in general. In the 
absence of the 0-term, because they are parity even they only create an even number of pion 
fields, and therefore will only contribute to valence quantities beyond NLO. The term with 
no pions from V ' s gives rise to a mass correction of baryons with at least one sea quark. 
These are also higher order than we are working. However, these operators will contribute to 
vertices in the one-loop graph contributing to the neutron EDM, replacing the tuq spurion 



insertions in Fig. 1 of Ref. [92]. 



The LO contribution to the neutron EDM comes from a loop diagram in \FT [93], with 



one interaction given by the LO HB^PT Lagrangian |94| , 95, 96, |97| , or Eqs. ( |2.30| ) and ( |3.6| ) 
in the MA theory. After performing the U{\)a rotation to remove the #-term from the QCD 
Lagrangian, the quark mass term becomes (for SU{A\2)) 

= diagimue^^de^^rrije^^m^^mue^^de^) , (3.7) 



16 For Wilson sea fermions, there will be similar operators but which only scale linearly in the lattice spacing. 



These are the generalizations of Eq. (3.3) to the mixed action theory. 
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with the operators in Eq. (3.6) picking up similar phases, and thus provide new contributions 
to the loop graphs for the neutron EDM. This is not the entire story however. The values 
of the phases, (j) u ,4>d etc., are determined by the vacuum stability of the pion potential, 
requiring there to be no single pion vertices, which in the partially quenched theory leads to 
the relations in the small angle limit |92]1 



m u (j) u = m d (j) d = rrijfij = mrfi , (3.8) 



— <t>d = —, . r , 4>j = 7 1 \ > <t>i = 7 • ( 3 - 9 ) 



m u {rrij + mi) ' rn d (mj + mi) ' (rrij + mi) ' (rrij + mi) 

At finite lattice spacing, there will be additional contributions to the sea-sea meson potential 
arising from chiral symmetry breaking operators which will also pick up phases under the 
above mentioned U(1)a rotation. These will be the same operators which provide additive 
mass corrections to the sea-sea mesons and their contributions to the neutron EDM (and 
vacuum stability) can be easily accommodated with additive corrections to the sea-quark 
masses in Eq. (|3.8|), which will depend upon the particular lattice action used in the sea- 
sector. For example, with Wilson fermions, the pions receive an 0(a) mass shift. Requiring 
there to be no single pion terms from the LO pion potential leads to 

m Z = m m + aW /2B , (3.10) 

where W is related to the chromo-magnetic condensate, which appears in Eq. ( |3,4| ) for ex- 
ample. This is defined in a similar fashion to the chiral condensate, B. For improved Wilson 
fermions, there will be a similar additive correction at 0{a 2 ) which can be determined from 
the meson potential given in Ref. p0| of a similar form, 



m 



iw 

3,1 



m m + a 2 W/2B . (3.11) 



The staggered sea quark masses are protected from multiplicative mass renormalization 
by the taste-5 chiral U(1)a symmetry of the staggered action. However, the staggered meson 
potential is not invariant under a taste-singlet U(V)a rotation and therefore there will be 
additive shifts to the vacuum stability condition proportional to the taste-Identity meson mass 
splitting, similar to the additive corrections to the topological susceptibility with staggered 



fermions [98|. This amounts to a correction of the sea quark mass of Eq. ( |3.8j ) for staggered 
fermions of, 

m f,i 9 ' = m m + a 2 Ar/2Bo . (3.12) 

Putting this all together, and specializing to the case of degenerate sea quark masses 
and lattice spacing corrections which begin at 0(a 2 ), the MA extrapolation formula for the 
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neutron EDM determined with our modified prescription and Ref. |)2| is then given by 17 
APQ) . 



d 



N 



eVm sea 
(4tt/ w )2 

e6a 2 



AF w ]n 




#J + 2 C(/X) 



1 



(3.13) 



In this equation, the sea quark mass, m sea is given by Eq. ( 3.11 ), ( |3.12| ), or the appropriate 
variation thereof for a given sea quark discretization method. Furthermore, it is simple to 
accommodate non-degenerate sea quark masses as in PQQCD p2[, in which case rh sea is really 



proportional to the product of sea quark masses included in the chiral Lagrangian, including 
the appropriate additive mass renormalizations discussed above. For the case of Wilson sea 
fermions, the only difference is that the second line of Eq. ( p.13 ) scales linearly in the lattice 
spacing, as opposed to the quadratic scaling given, and rh sea is given by Eq. (|3.10| ). These 
contributions arise from the mixed action generalization of Eq. (^^) . 

In Eq. ( [Tip , F n and Fj u are combinations of the coefficients of the operators in Eq. ( [2.30| ) 



and can be found in Eq. (30) of Ref. [92] and the new mixed action contribution has a 
coefficient 

Qj + qi ' 



F a 2 



9AOi a 



1 

3 + 
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Pa i \ ( Ota Pa \ 

T -(<&+«) {t + y) 



(3.14) 



The sea quark electromagnetic charges are given by qj and qi. We see that even with these 
new considerations, the form of this extrapolation formula is still independent of the type of 
sea fermions used, provided the leading lattice spacing effects are 0(a 2 ). There will be similar 
modifications to the extrapolation formula relevant to all quantities which are sensitive to the 
chiral anomaly. 



3.2 Twist-2 matrix elements of the nucleon 

Twist-2 matrix elements are related to moments of generalized parton distribution func- 



tions ||, [10 
elements in flO 
L06H and partially quenched 



101]. Chiral perturbation theory was first applied to forward twist-2 matrix 



103 , 104 ] then applied to off- forward matrix elements [105]. The quenched 
~~ 1[ versions followed subsequently. Our approach can be used 
to convert the leading meson, single nucleon and multiple-nucleon twist-2 matrix elements 
to their mixed-action versions. One of the more important twist-2 matrix elements is that 
related to the axial charge of the nucleon, which has recently been computed by the LHP 
Collaboration using a mixed action scheme ||. This is part of a more ambitious program to 
determine the structure of nucleons with lattice QCD and provides an important benchmark 
for the calculation of other twist-2 matrix elements (l^, [18|. The nucleon axial charge can be 



Using our arguments above regarding the vanishing of the neutron EDM if one of the sea-quark masses is 
zero, allows us to uniquely determine two of the counterterms in Ref. those being d = f = 0. Also note 
the opposite sign convention of our Lagrangian, Eq. (2.30) as compared to that used in Ref. [92[. 
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calculated with the nucleoli matrix element of the axial- vector current 



(3.15) 



which can be mapped into the heavy baryon chiral Lagrangian [94, |9l|. To determine the 
extrapolation formula for the calculation performed by LHPC, one can use the partially 
quenched formula worked out in Ref. [62], with the particular choice of extending the axial- 
charge matrix to 

r 3 = diag(l, -1,0, 0,1,-1). (3.16) 

Recently, the mixed action extrapolation formula for the neutron to proton axial matrix 
element was determined which in the isospin limit is equivalent to the proton-proton 
matrix element with the T3 current, and relevant for Ref. Q. In the mixed action EFT, the 
extrapolation formula is given at NLO by 



( N(p) I ( M ^ 5 I N(p) ) 



2 tip S n Up 



(MA) r , (MA) 

A pp T , 



PP 



(3.17) 



where 



(MA)j, 



pp — 9 A 



,91 



L{rh ju )- (l2g A + 2Ag\ + 16g 2 A9l + YIg A g\ + gf) 



L(m w )— (I6g\ + !7g A gi + gf) + 2g A (g A + gi) 2 A 2 ju -^—^L(m n ) 
^(4fiu + gi)g\ N K(m 7r , A,//) - ^g A - gi)g\ N K{m ju , A, fi) 
+ -//a \ ! 9 A + 7^9 AA ) J (m n , A, n) + 2g\ N (g A + ^9AA ) J {m ju , A, //) 



(3.18) 



and 



(MA), 



mlC m + A 2 u C^ + a 2 C a . 



(3.19) 



is given by local counter terms. Our formula is in agreement with that in Ref. [31], however as 
we will explain shortly, we have a slight disagreement with the analysis presented in Ref. |31] 
in the estimation of the size of the lattice spacing dependent corrections. Before discussing 
the relevance of this formula to the LHPC calculation, first we contrast this formula with the 
continuum ;\RT formula which was used to perform the chiral extrapolation in Ref. ll, 



Tpp — gA 



1 



2(g A + 2g 3 A )L(m n ) 



and 



+ -%^( 81 SA + 25gAA)J(m 7r ,A,fi) - ^-g A g 2 AN K(m 7T , A, fi) 



mlC m 



(3.20) 
(3.21) 
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In both the xPT and MA formulae, L, J and K are chiral logarithm functions defined in the 
literature, with 

L{m) =m 2 \n{^j , (3.22) 



and J and K can be found for example in Eqs (62) and (69) or Ref. [62] respectively, and 



their finite volume equivalents in Ref. [107]. There are two important distinctions between 
the MA and xPT formulae. The MA extrapolation formula has two more counterterms 
than the continuum formula, as seen by Eqs. ( 3.19 ) and ( |3.21| ), however with the tuning 
used by LHPC, Aj u ~ and therefore they require only one more counterterm. The lattice 
spacing dependent counterterm in Eq. ( [3.19| ) is required by scale invariance and also simply 
follows from the spurion analysis presented in Sec. 2.1 or from the arguments in Ref. [31]. 



Additionally, the loop corrections in the MA formula depend upon two different pion-nucleon 
axial couplings, gA and g\ as opposed to the one coupling, gA in the continuum xPT formula, 
Eqs. ([H|) and ( |3T20| ) . The MA formula also depends upon the mixed valence-sea mesons 



as well as the taste-identity staggered pion mass. However, the staggered taste splittings 
are well known [^] and the mixed meson mass renormalization has recently been calculated 
and is also known well |jl7f |. Therefore the MA extrapolation of the nucleon axial charge 
calculated by LHPC requires the determination of two additional unknown, unphysical terms 
as compared to the continuum xPT formula. 

Given the presence of these two unphysical LECs, and the limited amount of mass points 
in the LHPC calculation, one would like to estimate the size of the corrections to the contin- 
uum extrapolation formula due to the mixed action artifacts, to determine the impact these 
unphysical effects have in the extraction of the nucleon axial charge. This was taken up in 
Ref. pi] , and we do not repeat the analysis here, but we highlight a point of disagreement we 
have regarding the size of the lattice spacing dependent corrections. In Ref. [31], the lattice 



spacing counterterm, C a of Eq. (3.19), was varied in an uncorrelated fashion with the par- 



tially quenched pion-nucleon coupling, g\ and the mixed meson mass renormalization, which 
was not known at the time. This lead to a predicted error band of 0(200%) of the value of 
the nucleon axial charge, g\ as measured at m n ~ 350 MeV by LHPC (Table 1 of Ref. ||) 
which was almost entirely due to the mixed action lattice spacing dependent corrections, see 
Fig. 3 of Ref. pl|. This signals either a breakdown in the mixed action expansion for this 



observable or an overestimate of the errors in Ref. [31]. Given the quality of the lattice results 
and extrapolation performed by LHPC || and the small corrections the mixed action effects 
contribute to other quantities |25|, ^] we find the latter to be the more plausible explanation. 

At a fixed lattice spacing, by performing the chiral extrapolation there is no way to 
distinguish between the lattice spacing dependent counterterm, C a , the pion-mass indepen- 
dent contributions from ^ MA ^T pp and the LO contribution to the axial charge. Therefore, 
when LHPC performed their chiral extrapolation and determined a best fit value for chiral 
Lagrangian parameter gA, they were actually determining the linear combination 

~9A = 9A + a 2 C a (3.23) 
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where C a is a linear combination of the lattice spacing dependent counterterm, C a and loop 
contributions to ( MA 'Tpp which do not vanish in the chiral limit and are proportional to either 
a 2 A/ from the hairpin contribution or a 2 AMix from the mixed meson mass contributions. If 
we then assume a perturbative expansion (which breaks down for light enough pion masses, 
as the hairpin interaction diverges in the chiral limit |6^, |3T|), the uncertainty in this ex- 
tracted parameter is a much better estimate of the size of the lattice spacing artifacts in 
LHPC's calculation of the nucleon axial charge, and this is already contained in the error 
band presented in Figure 1 of Ref. [0]. In support of this estimation, it has been found with 
other quantities, that the general size of the mixed action artifacts with the LHPC mixed 
action scheme [|], ^] have been at the 1-5% level |29|. A full analysis of the mixed action 
corrections will involve the determination of the coupling g\ as well as making use of both the 
staggered taste-identity pion mass splitting [^5| and the recently determined mixed meson 
mass renormalization (l?]] which is beyond the scope of this work. We finally note that it 



was show in Ref. [31 that it is the uncertainty in the partially quenched parameter g\ which 



dominates the uncertainty in the mixed action corrections to the nucleon axial charge. 

More important to the LHPC program of calculating the structure of the nucleon 15, 18], 



one can use our prescription, the lattice spacing mass shifts calculated in Refs. |85|, 17] and the 



extensive use of the work of Detmold and Lin in Ref. [105], in which the partially quenched 
extrapolation formulae for all the forward twist-2 matrix elements has been determined to 
NLO in both finite and infinite volume, to determine the corresponding mixed action formu- 



lae. 18 Using the spurion analysis presented in Sec. 24., one can show that each twist-2 matrix 
element will have its own lattice spacing dependent counterterms which can be treated as 
a 2 C^ at NLO, as in Eq. ( 3.1 9| ) . As with the nucleon axial charge, the partially quenched 



pion-nucleon coupling g\ will likely play the largest role in the mixed action corrections be- 
cause it is precisely the linear combination (g^ + 51) 2 which is the coefficient of the hairpin 
interaction contributions to all the twist-2 matrix elements. Therefore it will be important 
to determine this coupling to have good control of the mixed action lattice artifacts. 

3.3 NN Scattering 

One of the greatest challenges facing the nuclear physics community is to determine the 
properties of nuclei from QCD. Given the success of lattice QCD with the meson and single 
nucleon sector calculations, it is natural to use lattice QCD to study nuclear systems. One 
of the complications however, is that nuclear physics is a finely tuned system. The deuteron 
binding energy for example B<i ~ 2.24 MeV, is much smaller than the scale set by pion 
physics. This makes it quite formidable to extract the deuteron binding energy in a lattice 
calculation out of the approximately 2 GeV rest mass of the proton and neutron. It turns out 
that by using heavier than physical light quark masses (corresponding to ~ 300 MeV), the 
two nucleon scattering lengths, which can be related to binding energy, tend to be a natural 



18 We note that the iso-vector twist-2 matrix elements were first determined in Refs. | |6l| , p2| ] for SU(6\3) and 
5(7(412) partially quenched theories respectively. 
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size and measurable with a box of L ~ 2.5 fm per side |7]]. 19 To further extract the physical 
scattering lengths, extrapolations to the continuum infinite volume limit and the physical 
quark masses are required. 

Recently, the NPLQCD Collaboration has determined the nucleon-nucleon scattering 
lengths in the 1 Sq and 3 Si channels using a MA lattice calculation with domain- wall valence 
quarks and the asqtad improved staggered MILC gauge configurations J5J. It is therefore of 
considerable interest to understand how the MA artifacts pollute the correlation functions. 
In Ref. [73], Beane and Savage have developed the partially quenched version of the two 



nucleon systems based on BBSvK power counting gig ]. In BBSvK, the channel follows 



the KSW power counting 71, 72| with one pion exchange entering as a perturbation at NLO 
while the 3 Si channel follows Weinberg's power counting [68, S9, 7C] with one pion exchange 



entering non-perturbatively and being resumed to the LO one pion exchange potential. The 
PQ effects (Aj u ^ 0) in the two- nucleon system arise from several important sources: the one 
pion exchange diagram also includes a hairpin interaction modifying the long-distance part 
of the potential [f73[|, the PQ corrections to the masses and couplings of the particles; qa, fw, 



m 7 



1 1 a \ 

and rriN, as well as two new PQ NN couplings, D^ B and D. 



)( 3 Si) 
J 2B ■ 



Using the prescription described above, we can easily incorporate the MA effects in two 
nucleon quantities once the corresponding PQ effects are known. Now we go through the list 
of PQ effects listed above and modify them to include the MA effects. The one pion exchange 
potential (OPE) is 20 



V^ P A E {r) 



<7l • V C?2 • V 



- - A 2 
2 r l ' r 2 / , n2 ^ju 

9A — {9A + 9l) 2^ 



(3.24) 



in the 1 Sq channel. The only modification from the partially quenched potential determined 
in Ref. |73| is Aj u — > Aj u . 21 The formula for the pion decay constant has been worked out in 
several places and we list it here for convenience, 



U = f 



1 



(4tt/)2 



In 



+ : 



777 



A 2 2 

{PQ), ^ia JMA\, . Ct 



V)-^+e M V) /2 



(3.25) 



19 There is some subtlety here. When determining the two-particle energy levels from the spatial volume 



110 



111 



dependence of the four-point correlation function, commonly referred to as Liischer's method 109 
|l!2| ], it is not the scattering length which determines the required size of the volume for the method to be 
applicable, but the effective range, generally set by the inverse pion mass in QCD. Therefore, even with 
unnaturally large scattering lengths as with the two-nucleon system at the physical quark masses, one can 
determine the infinite volume scattering parameters from the two particle interaction energy even when the 
scattering length is much larger than the finite spatial extent of the lattice |ll3fl . A smaller scattering length 
in the 3 Si channel is indicative of a larger binding energy for the deuteron making it easier to determine 
on the lattice. Also, the inverse box size, L _1 determines the splitting in energy levels of the two-particle 
system which is important for having well separated eigenstates. And lastly, one needs to make sure the 
exponentially suppressed volume modifications, which generically scale as e~ mxi are under control in the 



two-nucleon system [114], which is not as straight forward as with the two pion system [115]. 

20 The hairpin modification to Eq. (3.24) appears different to that in Refs. j?3|, [ii]], however this is simply a 
different convention for labeling the pion-nucleon couplings. Our convention is consistent with Ref. |62j]. 



Recall that A, 



^ju \a- 
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The MA formula for the nucleon mass can be found in Ref. J24J and gA, which we showed 
in the previous section, can also be found in Ref. |3l]]. The NN counterterms, following the 



spurion analysis presented in Sec. 2.1, should be replaced by 



D^ml + A 2 ju D^ + a 2 D^. 



(3.26) 



Making use of the work in Ref. [54|, we then find that in a MA theory with Ginsparg- Wilson 
valence fermions, the 1 S'o scattering length and effective range extrapolation formula are 



+ (ml - if 



A 2 (9A + gi) 2 M N 

JU 87T/I 



In 



+ 



47T 

1 



1? 



J_ 



(3.27) 



and 



M N 
2vr 



^- 1 fC 2 {ii) + 



9a M n ( 3 _ ^SL 
12vr/2 V mv 



+ 6- 



+ 



A j U (g A + gi) 2 M N 
ml 6vr/2 



2^-3^ 



where 7 is a /^-independent linear combination of \i and Cq (ji), the LO iV./V" interaction. 
Compared to the PQ case, the scattering length has one new lattice spacing dependent 
counter term at this order, D\ a (fi) while the effective range does not depend upon any new 
counterterms. The results in 3 S*i channel unfortunately do not have an analytic form because 
it requires solving the 3 Si- 3 Di coupled Schrodinger equation. However, the difference in 
the PQ and MA potentials follows the description given above, and can be determined from 
Ref. |)4|] with our prescription. 

Lastly, we comment that this discussion naturally extends to the hyperon-nuclear inter- 



(3.28) 



actions as well, for which the partially quenched theory has been developed in Ref. [117]. This 
is also very relevant as the first lattice study of the hyperon-nucleon interaction has recently 
been performed in a MA scheme as well [pi . 



4. Discussion 



In this work, we have proven that the new leading order meson operator allowed by mixed 
action theories, which is independent of the mixed lattice action, functions exactly as a 
meson mass operator for all vertices with (2A^ — 2) valence-valence mesons and 2 mixed 
valence-sea mesons. This proof, combined with the very nice features of mixed action effective 



field theories with chirally symmetric valence fermions 2E] has allowed us to construct a 
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prescription to determine many mixed action extrapolation formulae for valence quantities, 
through the leading loop order given the corresponding formulae determined in partially 
quenched chiral perturbation theory. Our prescription works immediately for quantities which 
do not depend on the term; we have further assumed that the mixed action theory has the 
same hairpin structure as the partially quenched theory, and that the valence quarks are chiral. 
In the case of a mixed action theory with more involved hairpin structure than the partially 
quenched theory, a more general prescription is undoubtedly possible. We have used the 
neutron EDM as an example which requires a modification of our prescription as it depends 
critically upon the 6 term, and have discussed how to modify our prescription for this quantity. 
With the recently measured mixed meson mass renormalization [17] and the well known 
staggered meson mass taste splittings [85], these mixed action extrapolation formulae can 
readably be applied to a host of physical quantities covering a broad range of hadronic physics: 
pion and kaon physics M, fill, M, M, ||, pi, M, M, E 



29|, baryon observables |61, [K 



120) , [I2ll , heavy meson observables |K]] , heavy hadron observables 



H, H> P arit y 

violation [123, 92 1, electromagnetic properties and transition matrix elements [124, 125, 126, 
127 ], structure functions |108| , 128 ], two nucleon systems [73], 64], hyper-nuclear systems [ 117 ] 
and constraints on beyond the standard model physics from hadronic contributions of Ab = 2 
and Ac = 2 observables 111 2£ 
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A. Chiral valence fermions on twisted mass sea fermions 

Here we present a few technical details of the mixed action chiral Lagrangian corresponding to 
chiral valence fermions and twisted mass sea fermions alluded to in the main text. We focus 
on the chiral theory for two degenerate valence and ghost flavors as well as a degenerate pair 
of twisted sea flavors. We will demonstrate that this theory behaves nicely in the following 
sense; the 0(a 2 ) mixed valence-sea fermion mass shift is given entirely by the CMi x -term at 
maximal twist; the hairpin structure of the theory is the same as in partially quenched xPT 
at maximal twist to the order we are working. These two facts follow naturally from the 
properties of graded algebras and the behavior of two flavor twisted mass QCD. We first 
summarize the discussion for twisted mass xPT which is relevant to our discussion. 
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The Euclidean Symanzik Lagrangian for twisted mass lattice QCD g| is (l3| 

C e ff = Cgi ue + V> lp+m + i^ 5 T 3 fi if) + bia^ia flu F flu ip , (A.l) 

where C is the gluon action, m = Z m (mo — m c )/a is the standard quark mass and [i = Z^fio/a 
is the twisted quark mass. From here one can construct the low energy chiral Lagrangian 
which is given at LO by |76], 77, 78] 



f 2 



where 



X = 2B (m + ir 3/ u) = m + ir 3 /i 



J —%t (d^d^) - J —ti ( X f S + £t X ) - ^tr (A+E + EtA) , 



and 



A = 2W a = a . 



(A.2) 



(A.3) 



The first observation to make is one can define a shifted spurion field such that the Lagrangian 
takes its continuum form, with 



X 



X + A = Me 



(A.4) 



where M' = w (m + a) 2 + ft 2 . The second observation is that because of the twisted mass 
term, the vacuum will no longer be aligned with the identity, which can be determined by 
minimizing the vacuum energy. One finds 



S = (0|E|0) = e^ 0T3 , 
such that the physical fields are given by an axial rotation from E, 

E = e m/2 Z phys e iul0Ts/2 . 



(A.5) 



(A.6) 



Expanding the LO and NLO Lagrangians about the vacuum, one then finds the automatic 
0(a) improvement of physical observables [^, [77], [78| . We now extend this analysis to the 
mixed action theory. Working in the isospin limit of the valence and sea sectors, the mixed 
action Lagrangian including the leading 0(a 2 ) operators is 



C = C str f^E^Et) - ^str (V E + £ f x 



2 C M i x str (T 3 £T 3 E tN ) - W' [str ( A f E + tfA) 



(A.7) 



where here the mass spurion includes the lattice spacing mass shift to the twisted sea mesons 
and is given by 



X 



t> -1-2x2 



o 2 



x2 



M'e iui0T3 
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al 



2x2 



(A.t 



2x2, 



J2x2, 
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with 



M' s = ^{m s +a} 2 +(i 2 , 

m v = 2B m va i, m s = 2B m sea , fi = 2B fj,, a = 2W a, (AS 



The singlet field is defined as 



/ 

$ = — In sdetS = sti(p ■ 
2i 



(A.10) 



Ultimately the singlet will be integrated out of the theory but it is convenient to keep around 
to determine the structure of the neutral propagators [|52]]. 

A. 0.1 Vacuum angle and meson masses 

We first address the shifted vacuum caused by the twisted mass term. It is straightforward 
to check that the vacuum energy is minimized by expanding about 

^2x2 

S = £o ^phys Co with Co = 



Ju T 3 /2 



(A.11) 



V 



1 



2x2, 



Expanding Eq. ( A.7 ) around this vacuum, one then finds the valence- valence, valence-sea 
and sea-sea pion masses are given at arbitrary twist and LO by (we have neglected terms of 
0(m q a) here which are proportional to cos a;) 



m, : 



1 



1 



m v + -M' + a A M ix + 



32W 



a cos uj, 



( mf±r = M' s + — a cos^ w, 



/ 



. 2 6W„2 2 

a sm uj H ^— a cos uj . 



(A.12) 



We then see at maximal twist, the valence-sea mesons only receive lattice spacing corrections 
from the a 2 AMi x term, in agreement with Eq. ( |2.23| ). 

A. 0.2 Hairpin interactions 

We must also address the hairpin interactions. In Ref. |23|] it was shown that in addition 



to the Lagrangian, Eq. (A.7), there are additional hairpin interactions which arise from the 
operators 

1st ^ 1 ( 



5C 



32 



str(V s (E-£t) 



(A.13) 



Expanding about the twisted vacuum, Eq. ( A. 11 ), this leads to an interaction 

5£ = -jssO 2 cos 2 uj str (Vs<f>) ■ 



(A.14) 
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Away from maximal twist, this interaction acts like a shift in the partial quenching parameters, 
Eq. ( 2.2SD , although the coefficient ~/ ss would need to be determined. However, we see that 
at maximal twist, this extra hairpin interaction is absent. In the notation of Ref. |pl| , we can 
derive the form of the flavor neutral propagators at arbitrary twist. First ignoring the extra 
hairpin, one can show the valence-valence neutral propagators are given by 



r , Ai ^'r ) ; 1 raid',- . , 



where 




Vval = -7T , (A.16) 



and 




l6W'a 2 ( 2 cos 2 w — sin 2 u sin 2 w i _ 

+ 72— -2 o 2 ■ 2 • ( A - 17 ) 

j z \ sin u) 2 cos oj — sin lo 



The valence-valence flavor neutral propagators, including the extra hairpin interactions of 
Eq. ( |Al4|) are then given at arbitrary twist by 



S VlV2 IP + l m 7r J + \—p r la -f ss J COS LO 



(A.18) 



* V1V2 p 2 + (m^ 1 ) 2 2 {p 2 + {ml lVl ) 2 ){p 2 + {ml 2V2 ) 2 ) 
Defining the fields 7r° = -^(ilu ~ Vd) an d V = T/jC^u + ^d)' one finds a * maximal twist 

r 1 r\ r - K±) 2 -("C) 2 /a 

^°- p 2 + (mr) 2' aIld yfi- b 2 + (mr) 2]2 • ^ i9 i 

It is interesting to note that it is not the mass of the neutral twisted mass pion which enters 
the numerator of the hairpin propagator, but rather a mass which is equivalent to the charged 
pion mass, which at maximal twist is free of even 0{a 2 ) corrections. 
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